On comparing sums of square roots of small integers * 
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Let k and n be positive integers, n > k. Define r(n, k) to be the minimum positive value 

of 
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I where ai, 02, • • • , a^, 61, 62, • • • , are positive integers no larger than n. It is an important 

, problem in computational geometry to determine a good upper bound of — log r(n, k). In 

this paper we prove an upper bound of 2'-'("/i°g") logn, which is better than the best known 
result 0(2^*^ log n) whenever n < ck log k for some constant c. In particular, our result 
' implies a subexponential algorithm to compare two sums of square roots of integers of size 

O ■ o{k\ogk). 
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In computational geometry, one often needs to compare lengths of two polygonal paths, whose 
nodes are on a integral lattice, and edges are measured according to the L2 norm. The problem 
can be reduced to the problem of comparing two sums of square roots of integers. Most work 
in computational geometry assumes a model of real-number machines, where one memory cell 
5^ \ can store one real number. It is assumed that an algebraic operation, taking a square root as well 
as a comparison between real numbers can be done in one operation. There is a straight-forward 
way to compare sums of square roots in real-number machines. But this model is not realistic, as 
shown in |6ii3. 

If we consider the problem in the model of Turing machine, then we need to design an algo- 
rithm to compare two sums of square roots of integers with low bit complexity. One approach 
would be approximating the sums by decimal numbers up to a certain precision, and then hope- 
fully we can learn which one is larger. Formally define r(n, k) to be be the minimum positive 
value of 

Ia/oi H h \Ai \/h 
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where oi, 02, ■ • • , a^, &i, &2, ■ ' ' > ^fc are positive integers no larger than n. The time complexity 
of the approximation approach depends directly on — logr(n, k), since an approximation of a 
sum of square roots of integers can be computed in time polynomial in the number of digits 
in the approximation. One would like to know if —logr(n, k) is bounded from above by a 
polynomial function in k and logn. If so, the approximate approach to compare two sums of 
square root of integers runs in polynomial time. Note that even if the lower bound of — log r(n, k) 
is exponential, it does not necessarily rule out a polynomial time algorithm. 

Although this problem was put forward during the 1980s |3|, progress has been scarce. In 
Q, it is proved that 

- log r(n, A;) = 0(2^'' log n) 

using the root separation method. This immediately gives us a polynomial time algorithm of 
comparing sum of square roots if k is fixed. Qian and Wang gave a constructive upper bound 

of r(n, k) at 0(n~^^+5), which corresponds to a lower bound 

— logr(n, k) = Q{k logn). 

They conjecture that — log r(n, k) = Q{n^~'^'' ^). 

There is a wide gap between the known upper bound and lower bound of — log r(n, k). Until 
the fundamental problem has been resolved, we can not even put the presumably easy problem 
such as Euclidean Minimum Spanning Tree problem in P, and the Euclidean Traveling Salesman 
problem in NP. 



1.1 Our contribution 

From the known upper bound of — log r(r;,, k), we conclude that there is a polynomial time algo- 
rithm to compare sum of square roots if k is fixed. In this note, we consider the case in the other 
end of the spectrum when k grows (almost) linear with n. 

Definition 1 An integer n is called square free, if there does not exist a prime p such that p^ 
divides n. 



It is well known that there are about 



7r2 



many square free integers less than n. If ai, 02, ■ " " ? o.k, bi, b2, - ■ ■ ,bk are distinct square free 
integers, then their square roots are linearly independent over the field of rational number Q. So 
it is possible that k and n are linearly related. This case is also practically interesting. We often 
need to compare paths whose nodes are on an / x / integral grid. The distance between the lattice 
points are square roots of integers of size 0{P). There are many nodes in the grid, and if we 
select a dense subset out of the grid points, we arrive in the situation where n is linear in k. 

We obtain a lower bound of difference of two sums of square roots. Our lower bound beats 
the root separation bound as long as n < ck log k for some constant c. The corresponding 
upper bound on — log r(n, k) becomes subexponential when n = o{k log k), or more generally. 
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if the square free parts of the numbers grow at rate o{k log k). Our bound implies a subexpontial 
algorithm, i.e., an algorithm with time complexity 2°('^Mogn, to compare two sums of square 
roots of small integers. The proof is also simple. 

We begin the presentation of our result by defining the notion of multiplicative generators. 

Definition 2 Given two set of positive integers A and B, we say that B multiplicatively generates 
A if any number in A can be written as a product of numbers from B with repetition allowed. 

It is easy to see that A multiplicatively generates itself, but for many sets, there exist much 
smaller sets which multiplicatively generate them. For example, all the square free number less 
than n are generated by the set of primes less than n, whose cardinality is 0(n/ logn). 

Theorem 1 (Main) Let Ci, C2, • • • ,Ck,di,d2,--- ,dk be positive integers. Let 

A = {ai,a2, - ■ ■ , flfc, ^1, &2, ■ ■ ■ , h} 

be the set of 2k positive square free integers. Assume that cfai < n for all 1 < i < k and 
d}bi < nfor all 1 < i < k. Let B be a set which multiplicatively generates A. Then 

IciVoiH \- Cky/a^ - di\/h dky/h\ > (2A;Vn)"^''''"'"\ 

Since A generates itself, so this result recovers the best known lower bound on r(n, k). In 
many cases, this result improves that bound, since \ B\ can be smaller than 1^41 = 2A;. It is possible 
that the cardinality of B can be as small as 0(log k), in which case, there is a polynomial time 

algorithm comparing two sums of square roots. 

Our result shows that the multiplicative structure of A affects the minimum possible value of 

Ici^/a^ + C2-s/a2 -\ 1- Cky/a^ — di^/b^ — d2\/% dk\/%\, which appears to be unknown 

before. In particular, we show that the root separation lower bound 2*^^'^^ logn of — logr(n, k) 
is not tight, at least, when n is linear in k. It is still possible that when n is much larger than 
k, the root separation bound becomes tight. Our result indicates that to achieve the root separa- 
tion bound, it is important to select the numbers Oi, 02, • • • , a^, 61, 62, • " " > ^fc such that they are 
pairwise relatively prime. 

2 The proof 

Let F — Q{xi,X2, • • • , Xm) be the function field over Q with indeterminate xi,X2, - ■ ■ , Xm- 
Consider a field extension K = F[yi, ^2, • • • , Vm]/ {Hi ~ ^i-, - ■ ■ ^Vm ~ ^m) of F. It is a linear 
space of dimension 2^" over Q(xi, ^2, • • • , Xm), one of whose bases is 

{Bs ^HvilSC {1,2,- ■■ ,m}}. 

ies 

The Galois group G of K over F has order 2"*. For any subset -S" of {1, 2, • • • , m}, define as & G 
recursively as follows: 
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1. If S — $, as is the identity element. 

2. If \S\ = 1, then 



-Vj if« = J 



3. lf\S\>l,as = U^esm■ 
We have (75' (^5) = (-Ij'^'^^'^s and G = {a^l-S C {1, • • -m}} 

Lemma 1 Le?{Q;5|5' C {1,2,3, • • • ,111}} be a set of 2"^ integers. The norm ofY^g^-^^ ^ ... m}^sBs, 
denoted by 

Nk/f{ ^ oisBs) 

5C{l,2,-,m} 

is a polynomial in Z[xi,X2, ■ ■ ■ , Xm]- 
Proof: By definition, 

nk/f{ yi "*^^) = E «^^^) 

5C{l,2,-,m} aeG 5C{1,2,- ,m} 

n ( E (^s<JsiBs)) (2) 

S'C{l,2,-,m} 5C{l,2,-,m} 

= n ( E (3) 

5'C{l,2,-,m} 5C{l,2,---,m} 

The norm must be an element in F = Q{xi,X2, ■ ■ ■ , Xm)- On the other hand, if we expand the 
product in the right hand side, it reduces to ^5c{i 2 ••• m} PsBs, where Ps ^ Z[a;i, X2,--- , a;^] 
for any 5" C {1, 2, • • • , m}. Hence we must have ^5 = for \S\ > 1. Thus 

Nk/f{ asBs) = /?0, 

SC{l,2,---,m} 

which is a polynomial in Z[a;i, X2, • • • , □ 
Define the polynomial 

SC{l,2,-,m} 

Now we are ready to prove the main theorem. 

Proof: (of the main theorem) Denote \B\ by m. Assume that B — {hi, h2, - ■ ■ , hm}- There 
is a natural ring homomorphism 
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by letting ip^i/i) = y/hl and 'ip{xi) = hi for 1 < i < m. 

Fix an order among all the subsets of {1, 2, ■ ■ ■ , m}. Define B'g = Ylies ^^'^ define as as 

{Cj if S is the first set such that aj = B'g 
—dj if S is the first set such that bj = B'g 
Otherwise 

We have 

ciVoiH hCfc^Ofc-c^iV^ dk^/h= ^ asB'g 

5C{1,2,--- ,m} 

and 

S'C{1,2,- - ,m} 5'C{1,2, -- ,m} 

because of the ring homomorphism. The integer 

since y/d^, y/d^, ■ ■ ■ , ^/d^, y/b^, y/b^, ■ ■ ■ , are linear independent over Q. So 

I n ( E (-i)'^'"^'"5B^)i>i 

5'C{1,2,--- ,m} 5C{l,2,---,m} 



Thus 



CiVoiH h CfcA/ofc - (iiV^i 4v&fc| (4) 



- n|5'|^0(E5C{l,2,....w(-l)"''"^'"^^^) 

□ 



The proof relies on the fact that the norm is a nonzero integer, thus has absolute value greater 
than 1. Every factor in the definition of the norm is not too large (less than 2kyfn in our case), 
so the smallest factor should not be too small. The technique has been used in several papers, for 
example, see The estimation depends primarily on the number of factors in the definition of 
the norm. 



3 A corollary from the main theorem 

Theorem 2 Let ci, C2, ■ ■ ■ , c^, di, d2,--- ,dk be positive integers. Let Oi, 02, ■ ■ ■ , a^, bi, 62, ■ ' ' ? ^fc 
be distinct square free positive integers less than m. Assume that cfai < nfor all 1 < i < k and 
dfbi < nfor all 1 < i < k. Then 

IciVOi H \- CfeVOfc - "iV^^i dk\/bk\ > [2k^/n) 
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Proof: It is well known that the number of primes less than m is 0{m/ logm). The set of 
primes less than m generates all the positive integers less than m. The theorem follows from the 
main theorem. □ 



Corollary 1 - logr(n, k) = 2°("/i°g") logn 

4 Conclusion remarks 

In this paper, we prove an upper bound of 2'^*^"/ log n for — log r(n, k), by exploring the fact 
that the algebraic degree of sum of 2k square free positive integers can be much less than 2^^. We 
suspect that 2*-'('=/^°sfc) logn type of upper bound holds for much large n, and leave it as an open 
problem. 
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